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Abstract 



A 2-level atom with degenerate ground state interacting with a quantum field is investi- 
gated. We show, that the field drives the state of the atom to a stationary state, which is 
^T} ! non-unique, but depends on the initial state of the system through some conserved quanti- 

ties. This non-uniqueness follows from the degeneracy of the ground state of the atom, and 
when the ground subspace is two-dimensional, the family of stationary states will depend on 
, a one-dimensional parameter. Only one of the stationary states in this family is a pure state, 

and this state coincides with the known population trapped state (zero population in the 
excited level \NC)). Another one stationary state corresponds to an equal weight mixture 
O ■ of the excited level 1 3) and of the coupled state | C) . 

m 

O . 

1 Introduction 

In the present paper we consider a 2-level atom with a degenerate ground state (or, equivalently, 



a 3-level atom with equal energies of the two lower levels). We prove that the interaction with 
radiation drives this atom to a family of stationary states, depending linearly on a one-dimensional 
parameter, which varies in explicitly determined interval. For a particular (extremal) value of this 
parameter the stationary state coincides with the coherent population trapped state, described in 

[l]-[5]- 

Our starting point are the papers [2], [4], which discuss coherent population trapping (CPT) 
in a 3-level A-system (i.e. a 3-level atom where only two transitions between the the lower levels 
1 and 2 and the higher level 3 are allowed, and the transition between 1 and 2 is forbidden). CPT 
is based on the preparation of atoms in a special coherent superposition of the two lower states. In 
[4] it was argued that the CPT process may be described in the basis of coupled and non-coupled 
states defined by 

|C) = ^(|1) + |2)) (1) 

|M7> = -^(|1>-|2» (2) 

where the two levels |1), |2), correspond to the hyperfine split 3Si/ 2 sodium ground state, coupled 
by the laser fields to a common excited level |3) within ?>P\/2 levels. 

In [2], [5] it is explained that this scheme, although highly simplified with respect to the real 
situation, nevertheless captures the main physical features of the phenomenon of CPT. 
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Coherent population trapping (CPT) consists in driving the atomic state to the superposition 
wave function \NC), often called population trapped or dark state cf. [l]-[5]. An atom driven 
in this state is transparent, i.e. it does not absorb the incoming laser radiation, in the sense 
that transitions from the trapped state to the excited |3) are forbidden [3]. The idea is that, at 
resonance and in the stationary state, all the atomic population is pumped into the non-coupled 
state, which is not excited by the laser radiation, so that the excited state population, hence the 
emitted fluorescent intensity reaches the minimum. 

In [6] trapped states were used to propose a scheme to utilize photons for ideal quantum 
transmission between atoms located at spatially separated nodes of a quantum network. 

Now, the process of driving the atom to the non coupled state is dynamical and therefore it is 
natural to try and apply the stochastic limit technique [7] in order to deduce a (non phenomeno- 
logical) master equation for which the \NC) state is an attractor. The present note essentially 
confirms this intuition, but also shows that the situation is more complex than this. In fact the 
\NC) state is indeed an attractor for the reduced dynamics, but not a global one: there exists 
a one dimensional parameter family of stationary states and the precise interval of this param- 
eter is determined in Theorem 4 below. This means that the set of atomic states is split into 
nonintersecting sets, one for each value of the parameter, and each of these sets is a domain of 
attraction for the state corresponding to the given value of the parameter. These stationary states 
generally are mixed. Only one of these states is pure and coincides with the known non-coupled, 
or population trapped state (2), investigated in [l]-[5]. 

This extends the effect of population trapping and predicts a dependence of the trapped state 
on the preparation procedure. This also extends our ability to control the quantum states of the 
atom. In fact, by preparing the initial state, switching the interaction with the field and (possibly) 
filtering the excited state, one can realize a switch between the \NC) and the \C) states. 

In the present note the role of the velocity of the atoms has not been investigated. It is natural 
to conjecture, from the analysis of [2], [5], that the selection parameter of the stationary states 
depends on this velocity. The explicit form of this dependence is now under investigation. Also 
the dependence on the initial state of the field is investigated. It is shown that the above scenario 
can be realized only in equilibrium or non equilibrium, but non vacuum, states. In the vacuum 
state new phenomena, such as quantum beats, may arise (cf. Remark 11 below). 

2 The master equation 

For the investigation of the dynamics of a 3-level A-system interacting with radiation we use 
the stochastic limit approach, [7]. In this approach one introduces a slow time scale t/A 2 , where 
A is a coupling constant for the interaction of the system with radiation. In the limit A — > 
the dynamics is given by Langevin and master equations, cf. [7], [8], which are unambiguously 
derived from the original Hamiltonian. For the mathematical discussion of Langevin and master 
equations see also [9], [10]. Evolution of the slow degrees of freedom of the filed in the stochastic 
limit approach was considered in [11]. 

We consider a 3-level system with degenerate (for example, hyperfine split) ground states |1), 
1 2) and the excited state |3). 

The interaction of the system with the radiation field is described by the Hamiltonian 

H = H S + H R + XHj (3) 
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where the system degrees of freedom are described by the Hamiltonian Hs'- 

ff s = e 1 |l)(l|+e 1 |2)<2|+e 3 |3)(3| 

where is the energy of the level \i) (note that S\ — e 2 ). 

The radiation degrees of freedom are described by the Hamiltonian 

H R = J2 I to(k)a*(k)a t (k)dk (4) 

i J 

where aj(fc) is a boson field with a mean zero gauge invariant Gaussian state characterized by the 
pair correlations 

{a*{k)a 3 {k')) = N t (k)5 l3 5(k - k!) (5) 

and i, j = 1,2 are the polarization indices. 

The interaction Hamiltonian Hi is defined as follows 

Hi = J Y,~9M^(k)D* a dk + h.c. (6) 

ia 

where a takes two values 1 and 2 and 

£>i = |l><3|, D 2 = \l)(2\ 

The free evolution of the interaction is equivalent to an effective free evolution of the boson 
field of the form 

e -«M*)-«) a .(jfe) 

where u = e 3 — S\ is the Bohr frequency, which is equal to the difference of energies of the two 
energy levels. 

By the stochastic golden rule [7] the rescaled free evolution of the field above, in the stochastic 
limit, becomes a quantum white noise b iu) (t, k), or master field satisfying the commutation relations 

[bUt, k), b) u ,{t\ k')] = 2n5^5 t ,5(t - t')5{uj{k) - uo)5{k - k') (7) 

and with the mean zero gauge invariant Gaussian state with correlations: 

(bt(t, k)b 3U ,(t\ k')) = 2n5^,5 i3 5(t - t')S(u;(k) - u)5(k - k')N t (k) (8) 

(M*, k)b* juj ,{t\ k')) = 2n5^5 l3 5(t - t')5(cu(k) - u)5{k - k')(N t (k) + 1) (9) 

The Schrodinger equation becomes a white noise Hamiltonian equation, cf. [7] , [8] which when 
put in normal order is equivalent to the quantum stochastic differential equation (QSDE) 

dU t = (-idH(t) - Gdt)U t ; t > (10) 

with initial condition Uq = 1 and where 

(i) h{t) is the white noise Hamiltonian and dH(t), called the martingale term, is the stochastic 
differential: 

/t+dt 
h(s)ds = J2 (D* a dB iaiV (t) + D a dB* aiv (t)) (11) 
iaui 
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driven by the quantum Brownian motions 

rt+dt r rt+dt 

dB ia0J (t) := dkg ia (k)b iuJ (T,k)dT =:/ b iu) (r, g ia )dr (12) 

Jt J Jt 

(ii) The operator G, called the drift, is given by 

G = E {(9ia\9ip)ZD* a Df3 + HS^) (13) 

where the explicit form of the constants (gialg^J, called the generalized susceptivities, is: 

r N(k) + 1 

(9ia\9ip)Z = ~i J dkg ia (k)g i(3 (k) ' _ - - — (14) 



= nj dkg ia (k)g l( i(k)(N t (k) + l)5(Lu(k)-Lu)-tP.P. j dk g ia (k)g i0 (k) ^ ^ 

r N(k) 

(9ia\gip)u =- i J dk 9ia(k)gip(k) ^^ -j— - - (15) 

= 7T y dk g ia (k)g if 3{k)Ni{k)5{uj{k) - u) - iP.P. ^ dk g ia (k)g^(k)-^ 



Njjk) 
' uj{k) — UJ 

We will use the notations 

Re (ft a |^)+ = 7r y dkg ia (k)g il3 (k)Ni(k)5(uj(k) - u) (16) 

Im (^|^)+ = P.P. / dkg~fk)g ip {k) ^ (17) 

(9ia\9ip)t = Re {9ia\9ip)Z + ilm {9i*\9ip)Z ( 18 ) 
{gia\gip)i = Re (fi'ialfi'i/?)^ - ilm (gia\gip)u (19) 

Note that for a = (3 the values (16), (17) coincides with real and imaginary part of the 
generalized susceptivities, but the are and not necessarily equal to real and imaginary parts for 
a ^ (3. The values in (18), (19) are related in the following way: the complex conjugation of 
(gia\gip)u is equal to (^|^ ia )+. 

We also use the notation 

(9 a \9f})u = J2(9ia\9ip)u (20) 



Remark 1. Typically the (?j a are matrix elements (cf. the description in section 4.9.3 of [7]). 
Therefore their dependence on the index a is often unavoidable. Therefore we will develop the 
theory, as far as possible, keeping this dependence explicit. In some cases, e.g. some particular 
classes of 3-level atoms, the assumption that the formfactors g,i a do not depend on the index a, is 
justified. In this case the formulae simplify and are easier to interpret. Thus situation is described 
in the Section 3 below. 
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Remark 2. Note that if the expectation of number operators Ni(k) in the reference state 
depends only on the dispersion uj{k) (in this case we will denote this value N(uj)), then we have 
the identity 

Re(g\g)+ N(u>) 
W Re(g\g)- N(u) + 1 W 

which shows that this quotient of generalized susceptivities does not depend on the formfactor g 
and is a natural non equilibrium generalization of the Einstein emission-absorption coefficient (cf. 
section 5.9 of [7]). 

The master equation for the reduced density matrix p(t) of the system for a general discrete 
system with the dipole interaction in the stochastic limit approach was found in [8]. For the 
considered degenerate 3-level A-system it takes the form 

^ = Im ( 9jl \g n Up, |3)(3|] - i Im { 9jl \ 9 ^[p, 

+2Re( gjl \ 9jl )Z (p 33 |l)(l| - \{p, |3)(3|}) + 2Re( 9jl \ 9jl )+ (p n |3)(3| - ±{p, |1)(1|} 
+ (i Im (gMZlP, |3)(3|] - i Im ( 9j2 \ 9j2 )+\p, |2)(2|] + 



f | -i Im (j& 2 |<&i)>, |2)(l|]+2Re(^ 2 | ftl )JP33|l)(2| + 2Re(^ 2 | ffjl ): (pi 2 |3>(3| - Up, |2)(1|} 



+2Re(g j2 \g j2 )Z [p 33 \2)(2\ - -{p, |3)(3|}j + 2Re ( 9j2 \ 9j2 )t ^ 22 |3)(3| - -{p, |2)(2|} 
+ Im ( 9jl \g j2 )+[p, |1)(2|] + 2Re( 9jl \ 9j2 )ZP3s\2)(l\ + 2Re ( 9jl \ 9j2 )+ fp 21 |3)(3| - ±{p, |1)(2|} 

2' 

(22) 

where, as usual [a, 6] = a& — ha and {a, 6} = ab + 6a. 

One of our main results is the following separation of the density matrix into parts correspond- 
ing to invariant subspaces of the evolution. 

Lemma 1. The vector space H(3) of the Hermitian 3x3 matrices is the direct sum of two 
subspaces, V , V\, which are invariant under the evolution, defined by (22): 

if (3) = V ®V 1 

A linear basis ofV is given by {|2)(3|, |3)(2|, |3)(1|, |1)(3|}. Any matrix in this space decays 
exponentially to zero under the reduced evolution if the real parts of generalized susceptivities 
(16) (for the indices a = 1, 2 and (3 = 3 and vice versa) are non zero. 

A linear basis ofV 1 is given by {|2)(1|, |1)(2|, |3)(3|, |1)(1|, |2)(2|}. This space contains all the 
stationary states for the evolution. 

Proof Direct verification from the right hand side of (master3). 

Remark 3. Notice that the space V\ = C|3)(3| © M, where M is the 2x2 matrix algebra 
generated by |1)(2|, is itself a *-algebra. 
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Remark 4. From Lemma 1 we deduce that the evolution of the density matrix p(t) can be 
split into the sum of two evolutions po(t) and pi(t), where po(t) is an off diagonal matrix and pi(t) 
is a density matrix. 

/ p 32 (t) p 33 (t)\ fp 33 (t) 
p(t) = p {t) + Pl {t) = \p 23 (t) + p 22 (t) p 21 (t) 

\pi3(t) J V p 12 (t) p n (t) 

Moreover, ||po(^)|| < e ~ ct , where 2c = min Re (gj a \9ja)ti J ' = 1, 2, a = 1,2, i.e. the off-diagonal 
part of p(t) (in Vo) decays exponentially whenever c > 0. 

3 Stationary states for the master equation 

In the present section we will describe the set of stationary states for the evolution, generated by 
the master equation (22). By lemma 1 the invariant states of (22) belong to space V\. On this 
subspace equation (22) reduces to the following system of three differential equations 

dp 22 {t) 



dt 

dpnjt) 
dt 



= 2Re {g 2 \g 2 ) UJ P33 - 2Re (g 2 \g 2 )^p 22 - (gi\g 2 ) ul p 2 i - (gM^Pu (23) 



= 2Re (gMuPw - 2Re (gM^pu - (g 2 \gi)+p 2 i - {g 2 \gi) u) pi2 (24) 



dp 12 (t) 



9i\9i)u + (^2^2)^)^12 - (9i\92) u Pu - (92\gi)ip22 + 2Re(g 2 \gi) LU p 33 ( 25) 



dt 

which together with the normalization condition 

P11 + P22 + P33 = 1 

the conjugation rule 

P*12=P21, Pll,P22,P33 £ R 

and the conditions of positivity of the density matrix discussed in the following Lemma, form the 
set of equations determining the evolution of density matrix. 

Lemma 2. The Hermit ian matrix 

/p 33 \ 
P = P22 P21 

V P12 pii) 

is a density matrix iff the diagonal elements satisfy 

Pll + P22 + P33 = 1, Pll,P22,P33 > (26) 

and the off-diagonal elements satisfy 

P12 = P21, IP12I 2 < P11P22 (27) 
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Lemma 3. When the susceptivities (g^ga^, do not depend on a, (3 (in this case we denote 
them (glgju) the system (23)-(25) of linear equations determining the evolution of the atom has 
the conservation law 

Pu(t) + p 22 (t) = p 12 {t) + p 21 (t) + C; V* (28) 

Moreover, if Re(g\g)~ > 0, then p\ 2 {t) + Pziif) converges exponentially in time to the stationary 
value 

, = me{g\g)-{l-C)-2Re{g\g)+C = 1-C-CRJ2 
912 P21 ~ A{Re{g\g)- + Re{g\g)l) 1 + R w 1 J 

where C is a real constant. 

Proof The conservation law (28) follows from the identity 

dpujt) dp 22 (t) = dp l2 (t) dp 21 {t) 
dt dt dt dt 

It implies that for fixed C the evolution of the system is characterized by the real function of time 
P\2{t) + p2i{t) which we denote 2s(t): 

s(t) = l(p 12 (t)+p 21 (t)) 

With this notation the system (23)-(25) implies that 
ds(t) 



dt 



-4 (Re (g\g)- + Re (g\g)+) s(t) + 2Re (g\g)-(l - C) - Re {g\g)+C (30) 



If Re(g\g) UJ > 0, then equation (30) implies the exponential decay of s(t) to the stationary value 
(29) and this proves the lemma. 

Remark 5. Note that the condition Re(g\g)~ > means that 



g{k)g(k)5(u(k)-u)dk ^ 

which is automatically satisfied when the support of the formfactor g{k) intersects with the reso- 
nant surface cu(k) = uj. 

The stationary solution of the system (23)-(25) is determined by the system of equations 



(g 2 \g 2 )^(p u + p 22 ) + 2Re (g 2 \g 2 ) u p 22 + (gM^p*^ + (gi\g 2 )+Pi 2 = 2Re(g 2 \g 2 ) UJ (31) 
2Re ( gi \ gi )Z(pii + P22) + 2Re (gM+pu + (<? 2 |<7i)+Pr2 + G? 2 |<?i)>i2 = 2Re ( gi \ gi )Z (32) 

((gi\9iU + P12 = -(gMuPu - (g2\gi)ip22 + 2Re (g 2 \gi)ZP33 (33) 

Remark 6. For different formfactors g a (k) the system (23)-(25) may have different behaviors. 
In the generic case for g\ ^ g 2 the stationary solution is unique. For instance when g\ is orthogonal 
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to g 2 (in the sense of the bilinear form (<?i|<? 2 )J), then the determinant of the system (31), (32) 
reduces to 

- (2Re (<? 2 | <? 2 )+2Re ( 9l \ 9l )~ + 2Re (g 2 \g 2 ) ~2Re ( 9l \ 9l )+ + 2Re (g 2 | (? 2 )+2Re (<?i | <?i)+) 
and whenever this determinant is non-zero, the solution is unique. 
When g 1 = g 2 the solution is non-unique due to Lemma 3. 

Now we are ready to formulate the following theorem describing the structure of the stationary 
density matrices. 

Remark 7. If for a, (3 — 1, 2 

{9a\gp)Z = 

in particular, in the Fock case, the stationary solutions of (23)-(25) (neglecting the trivial case 
when also (g a \g/3)Z = 0) i s characterized by the single condition 

Pll + P22 = 1, Pll,P22 > 

so that P33 = and p\ 2 is arbitrary and subject only to the constraints (27). 

Theorem 4. For {g a \gp)^ not depending on a, (3 and when 

Re(g\g)+>0 (34) 

the system (31 )-(33) of linear equations determining the stationary state of the atom possesses a 
family of solutions parameterized by the one-dimensional parameter: 




P= p 9 s I (35) 
where, in the notations (21 ) 

2fle(g|g) + (l + 2a) (l+jgjg, 

Pe 4Re(g\g)- + 2Re(g\g)+~ 2 + R„ { ' 

= 2Re{g\g)z-2Re{g\g)ts = 1 - sR^ 

99 4Re(g\g)- + 2Re(g\g)+ 2 + R„ 1 } 
The admissible values of the parameter s are precisely those for which 

1 1 / Reiglg^Y 1 1 , N 

= 1+ 0. - >*>"o ( 38 ) 



2(1 + ^) 2 V Re(g\g)-J ~ ~ 2 

Moreover, if (34) is satisfied, the solution of the system (23)-(25) converges, ast^oo, to the 
stationary state (35). 

Proof If g± = g 2 = g, then (31)-(33) take respectively the form: 

2Re (g\g)Z{pn + P22) + 2Re (^|^)+p 2 2 + (<?|<?)>2i + (g\g)+Pi2 = 2Re (g\g)Z (39) 
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2Re (g\g)Z(pii + P22) + 2Re (g\g)+ Pll + (g\g)+ P 2i + (g\gfipu = 2Re (g\g)Z (40) 



2Re {g\g)Zpi2 = -(g\g)ZPu ~ {d\g)tp22 + 2Re {g\g)ZPw (41) 
Taking the differences of (39), (40) and of (41) and its conjugate, we obtain 

2 (Re (g\g)+) (P22 - P11) + 2t (im (g\g)+) (p 21 - p 12 ) = 

2 (Re (g\g)+) (p 21 - p 12 ) + 2* (im (g\g)+) (p 22 - p n ) = 

In the following we will not indicates the brackets at Re(g\g) and lm(g\g). 
Taking the sum of two equations above and dividing by two, we get 

{g\g)^{p22 - pn + P21 - P12) = 

If (g\g)u 7^ 0, then since p 22 — pn is real, and pi 2 — p 2 i is imaginary, we obtain 

P22 = Pn, P21 = P12 (42) 

Then, the sum of (39) and (40) takes the form 

2Re (g\g)+(p n + P22 + P12 + P21) = 4Re {g\g)ZPw (43) 

Equations (42) and (43) imply that any stationary density matrix must satisfy the following 
condition: 

2Re (g\g)+(pn + P12) = 2Re (^|^)jp33 (44) 

Since under the condition (42) the equations (39) and (40) coincide, equations (42), (44) describe 
the general stationary solution for (23)-(25). Using (42), (44) and (26), we obtain 



2Re(g\g)Z-2Re(g\g)+p 12 2Re (g\g)+ + 4Re (^)+Pi 2 

4Re(g\g)- + 2Re (g\g)+ ' P33 4Re (g\g)- + 2Re (g\g)+ 



In particular, p 12 must be a real number. From (45) and (27) one sees that the positivity of the 
density matrix is equivalent to inequalities 

1 / ReMoV+V 1 1 

- 1 + — ^11^. > p > — 46 

2 I Re(g\g)-J ~ Pl2 ~ 2 K ] 

Conversely, taking any real value of pi 2 satisfying (46) and determining p n and p 33 by (45), one 
obtains a stationary state for the master equation (23)-(25). 

Let us now prove the convergence of the system to a stationary state. The system (23)-(25) 
implies 

j t {p22 - P11) = -2Re (^|^)+(p 22 - pn) + 2% Im (g\g)+(p 12 - P21) (47) 

^(Pi2 - P21) = -2Re {g\g)t{pi2 - P21) + 2% Im (#|#)+(p 22 - pn) (48) 
Adding these two equations we see that 

P22 - P11 + P12 - P21 = const e<- 2Re ^ +2i Im toWi) (49) 
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For the case ~Re(g\g)^ > the linear combination (49) converges exponentially to zero. Since 
P22 — Pu is real and p\ 2 — P21 is imaginary, we obtain that (49) converges to the state where 
P22 = P11 and P12 = P21 (and therefore real). 

Then, applying Lemma 3, we get that p\ 2 = p 2 i, P22 = P11 and p 33 converge to stationary 
values, which are controlled by the stationary value s = \{p\2 + P2i)- 

This finishes the proof of the theorem. 

Remark 8. Note that if Re(g|g)+ = and lm.(g\g)+ ^ then (49) implies that the system 
does not converge to a stationary state but has an oscillatory behavior. 

Since the generalized susceptivities are given by the expression 

Re(^| ffi )+ = nj \gi(k)\ 2 Ni(k)5(uj(k) - u)dk 

Re (3^) J =nf Mk)\ 2 (N t (k) + l)6(u(h) - u)dk 
Re(</|#)± = J2 Re (9i\9i)» 

i 

It follows that one has inequality 

Re(g\g)- >Re(g\g)+ 

One can see that for high intensity of radiation, i.e. when N^k) >> 1, one can put = 1- 

In this case the solution (35), (38) will be simplified as follows 



/if* \ 
i=s- s 

\ s ±=v 



1 1 

- > s > — 
4 - - 2 



The most interesting states correspond to the extremal values of the parameter p 12 - The 
minimal value of p i2 is — |, which correspond to the density matrix for the pure state \NC): 

Pmin = S I 1 -1 I = \NC)(NC\ = \ (|1)(1| + |2)(2| - |1)(2| - |2)(1|) 




where the vector 



|M7> = -^(|1>-|2» 



is exactly the coherent population trapped state (2), discussed in the literature [l]-[4]. In the 
same approximation the maximal value p 12 = \ corresponds to the density matrix 



Pmax 



x /2 0. , 

7 1 1 =^|3)(3| + -|C)(C| 

yo 1 1 - : 
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This state is mixed, but the state of the reduced system corresponding to levels |1) and |2) is pure 
with the state vector 

|C) = i=(|l> + |2» 

which coincides with the coupled state (1). 

Thus the application of the stochastic limit approach allows to generalize the coherent pop- 
ulation trapping phenomenon. We see that the family of stationary density matrices realizes a 
continuous interpolation between the coupled and the non-coupled state. 

Remark 9. To distinguish experimentally different population trapped stationary states, one 
can measure the following observable 

A=|1)(2| + |2){1| 

which for instance may describe the interaction of the hyperfine split levels with a magnetic field. 
In fact the measurement of A in the generic stationary state gives 

tr pA = + P21 = 2s 

and the different trapped states give different mean values of A. 

Remark 10. Equations (37), (38) imply that the minimum ground state population is achieved 
when the parameter s is maximum, i.e. 

ruin = 1 ( 1 _ Ru, \ = 1 1 = 1 N(U) + 1 1 

Pg 2 + RuK 2 + 2i?J 2 1 + R w 2 2N{uo) + 1 4 

Since the ground level population is 2p g , it follows that in any stationary state at least 1/2 of 
the population is in the ground level. On the other hand the above chain of identities shows, in 
that the region of high radiation intensity N(u) » 1, the estimate p g nm = | is almost exact 
and therefore we can conclude that, in this region, for any stationary state, the population of the 
excited state is about 1/2. 

Remark 11. Consider the regime when there is no decay to the stationary state and we have 
the oscillations. This regime is possible when Re (gj a \9j/3)t = 0- We consider again the case when 
the susceptivities {g a \gp)^ do not depend on a, j3 and all Re {gj a \9jp)u — 0- 

By Remark 5 after Lemma 3 it is natural to assume that Re(g\g)~ > and there is a con- 
vergence of s(t) = \ (puit) + P2i(t)) to its stationary value — | < s < \. Analyzing the system 
equations for the density matrix, one can check that in the considered case the dynamics in the 
invariant subspace V\ is described by equation (49), which takes the form 

p 22 - P11 + P12 - P21 = const e 2U Im ^ (50) 

where pij are complex numbers satisfying Lemma 2. 

This kind of pure oscillatory behavior without damping is related to the quantum beating. 

When Re (g\g)Z > 0, the off-diagonal matrix elements P13, P23 decay exponentially by Remark 
4, cf. [8]. We see that in the regime Re (g\g)Z > 0, Re {g\g)t = (which is satisfied, for instance 
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in the Fock (vacuum) state) the behavior of the 3-level degenerate A-system for large times is 
described by the oscillations (50), when p i3 = p 23 = p 3i = p 32 = 0ands(t) = i(pi 2 (t)+p 21 (t)) = 
const, — \ < s < \. 

In conclusion: in the present paper we investigated the interaction of an atom with a degenerate 
ground state with a quantum field. We find (under natural conditions for the formfactors) , that the 
evolution drives the atom exponentially to a stationary state. This stationary state is not unique, 
and the family of stationary states may be parameterized by a one-dimensional parameter. For 
a special (minimal) value of this parameter the obtained stationary state is pure and coincides 
with the population trapped state, known in the literature [l]-[4]. The obtained results show the 
possibility of emergence of mixed stationary states, which continuously interpolate between the 
coupled and the non-coupled states. This difference can be experimentally detected. 

In the case of special states (the Fock state) also the oscillatory behavior (50) is possible. 
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